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Abstract 

The in-plane optical conductivity has been successfully employed to obtain 

information about the coupling of the charge carriers to the spin degrees of 

freedom in the high T c oxides. We investigate how this inelastic scattering 

affects the out of plane charge dynamics. We consider both coherent (in-plane 

momentum is conserved) and incoherent (no momentum conservation) c-axis 

transfers and find that the two cases give quite a distinct c-axis conductivity 

as a function of energy u. Comparison of our theoretical calculations with 

the available data does not allow a definitive conclusion, but a momentum 

dependent coherent matrix element characteristic of the Cu02 chemistry is 

favored with the possibility of a subdominant incoherent contribution. 
74.20.Mn 74.25.Gz 74.72.-h 
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I. INTRODUCTION 

Carbotte et a/.i recently suggested that a signature of the 41 meV spin resonance which 
is measured in spin polarized inelastic neutron scattering experiments^ in the superconduct- 
ing state of optimally doped YBa2Cus06.95 (YBCO), is also seen in the in-plane infrared 
optical spectrum as a function of energy uo. Following this initial suggestion Schachinger 
and Carbotteifl proceeded to analyze optical spectra in other high-T c compounds and found 
that similar spin resonances develop in many of the cuprates but not in all. These authors 
provided a simple technique0iHl to extract from infrared data an approximate electronic 
carrier-spin fluctuation spectral density which not only gives a picture of the position and 
width of the spin resonance, but also gives an absolute measure of its coupling to the charge 
degrees of freedom. More recently, consideration of in-plane optical dataB in YBCO at several 
temperatures reproduced well the temperature evolution of the 41 meV resonance! measured 
in neutron scattering experiments and gave a carrier-spin fluctuation spectral density which 
showed considerable variation with temperature. It consists of a broad spin fluctuation 
background which extends, in frequency, to several hundred millivolts and which persists 
in the normal state. Superimposed on this background is a spin resonance contribution 
which exists only in the superconducting state for optimally doped YBCO and grows as T is 
lowered. The growth of the resonance reflects modifications to the spin spectrum connected 
with the development of the superconducting state as the temperature is lowered below T c 
(the superconducting critical temperature). It is found that feedback effects increase the 
stability of superconductivity and, within an Eliashberg formalism for the superconduct- 
ing state, lead directly to a ratio of the d-wave gap amplitude to the value of T c which is 
considerably larger! than the BCS value of 4.3 in agreement with experiment. 

Other observed properties of the superconducting state can also be understood directly 
within the same framework. Most prominent among these are (1) agreement with experiment 
for the value of the zero temperature condensation energy per copper atomjl^lill (2) the fact 
that only about one third of the total optical spectral weight condenses into the superfluid 
density at zero temperature.lll This arises because it is mainly the coherent part of the charge 



2 



carrier spectral density which condenses while the incoherent boson assisted background is 
largely unaffected by the condensation. (3) The large peak observed in the microwave 
conductivity around 35 K0 is understood as due to the collapse of the electronic scattering 
rates0 which results from a low frequency gaping of the spin fluctuation spectrum due 
to the onset of superconductivity. (4) There is also a corresponding peak in the thermal 
conductivity which the calculations explain.^! 

In this paper we wish to extend the work to the c-axis optical conduct ivitySHll to see if it 
too can be satisfactorily understood within a generalized Eliashberg formalism with kernels 
determined from in-plane optical conductivity data.ll! The out-of-plane optical response has 
been the subject of considerable recent interestUHl due in part to the observation of an 
important violation of the conventional Ferrell-Grover-Tinkham^l (FGT) sum rule. This 
sum rule is obeyed by conventional superconductors but not by some high T c cuprates in the 
direction perpendicular to the CuO"2 planes. In this case the missing spectral weight under 
the real part of the optical conductivity which disappears on entering the superconducting 
state, does not necessarily equal the superfluid density which condenses into a delta function 
at zero frequency u = in the real part of the optical conductivity. In this paper we will 
concentrate on the frequency dependence of the c-axis optical conductivity a c (T,u) at low 
temperatures T in the superconducting state for uj > 0. 



II. FORMALISM 

The out-of-plane conductivity a c (T,u) at temperature T and frequency uj is related 
to the current-current correlation function II C (T, iv n ) at the boson Matsubara frequency 
v n = 2nirT, n = 0, ±1, ±2, . . ., analytically continued to real frequency u, and to the c-axis 
kinetic energy 

(i) 

with e the charge on the electron and d the distance between planes in the c-direction. In 
terms of the in-plane thermodynamic Green's function G(k, iu n ) and for coherent hopping 



<?c(T,u) = - 

! , J 



n c (T, 



u + iO + )-e 2 d 2 (H c ) 
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t_i_(k) perpendicular to the CuC>2 planes 

U c (T,w n ) =2(ed) 2 TY / T, t l( k ) tT {^G(k,tu; m )f G(k, i u; rn + iu n )} (2a) 

and 

(H e ) = 2r^^ti(k)tr{f 3 G(k,ia; m )f3G(k,ia; m )} (2b) 

UJm k 

In Eqs. (0) the 2x2 Nambu Green's function G(k, iu m ) describes the in-plane dynamics 
of the charge carriers with momentum k in the two dimensional Cu02 plane Brillouin zone, 
iu n = i{2n + 1)ttT for temperature T, n — 0, ±1, ±2, . . ., and is given by 

/Vi ■ \ tu>(iun)f + Ckr 3 + A k (^ n )fi 

G(k, zw„) = — - 2 — - , (3) 

where the f's are the Pauli 2x2 matrices, Ck is the band energy of the charge carriers as a 
function of their momentum k, A^(iuj n ) is the renormalized gap and u(iu n ) the renormalized 
Matsubara frequencies. In our model these quantities are determined as solutions of Eliash- 
berg equations suitably generalized to describe a ci-wave superconductor in terms of a charge 
carrier-spin fluctuation spectral density I 2 x{^) which we have determined previously from 
consideration of the in-plane optical conductivity.BUfl The methodBIHzl is to use the second 
derivative ^^[ UJT ^ 1 { UJ )} °f the optical scattering rate r _1 (u;) in the superconducting state 
as a first estimate in the construction of a model for the underlying charge carrier - spin 
fluctuation spectral density 7 2 x(w). 

In Eqs. (0) the out-of-plane matrix element t±(k) can depend on the in-plane momentum 
k. Models have been summarized in the recent preprint of Sandeman and Schofield0 who 
refer to previous literatureJ^H^ A possible choice is £_i_(k) = t±, a constant. But, consid- 
eration of the chemistry of the Cu02 plane and of the overlap of one plane with the next, 
suggests a form £j_(k) = cos 2 (20) where <fi is the angle of k in the two dimensional Cu02 
Brillouin zone for the plane motion. This matrix element eliminates the contribution from 
nodal quasiparticles entirely from the c-axis motion. 

For incoherent impurity induced c-axis charge transfer Eqs. @) are to be modified. After 
an impurity configuration average we obtain 



II C (T, iu n ) = 2(ed) 2 T ]T E { f °^ k > ^™) f oG(k, zu; m + %v n ) } (4a) 

m k,k' 

(H c ) = 2TEE^{ f 3G(k,zo; w )f 3 G(k,2a; m )}, (4b) 



m k,k' 



with V^ k / the average of the square of the impurity potential. If the impurity potential was 
taken to conserve momentum, which it does not, we would recover Eqs. (|2|). Various models 
could be taken for V£ k ,. Here we use a form introduced by Kimll and Hirschfeld et a/Ji-3 



Vk,k' = \Vo\ 2 + \Vi\ 2 cos(20) cos(2<//) (5) 

with <p and <$' the directions of k and k' respectively. 

As written Eqs. (§) and (§) involve Green's functions in the imaginary frequency Mat- 
subara representation. To obtain the conductivity <J C (T, u) as a function of real frequency an 
analytic continuation in Eq. flU) to real frequencies is needed. This could be done by Pade 
approximates, but here we analytically continue the entire equations and work with real 
frequency Eliashberg equations. For a <i-wave superconductor with gap Ak = Ao cos(2</>) 
on a cylindrical Fermi surface where cf) is the polar angle in the two dimensional Cu02 
Brillouin zone, the basic equations for A and Q which include inelastic scattering due to 
a boson exchange mechanism between the charge carriers and described by the spectral 
density J 2 x(c<j)ill! take on the form 

Af-L-A^ • T (o^Wf ■ m/ A^ m ;0Qcos(20Q \' 

A(z/ + to; (p) = trcl g cos(20) |A(z/ — iuj m ) + \{y + ioJ m )\ ( = > 

m=o \ y u) 2 {iuo m ) + A 2 (iu m ; <fi') I 

+ra Idz cos(2^ W H z) + /(*-*)]( A( y -, + ^)cc(2^) \' 
4o \ v / cu 2 (z/-2 + i5) + A 2 (z/-^ + i5;0 / )/ 

and in the renormalization channel 



V + i A ) = v + mT V] [A(z/ — iu m ) — \{u + ioj. 



U0[lLO ri 



OC 



m=o \yuj 2 (iuj m ) + A 2 (zu; m ; </>')/ 

■i- I dzI 2 X (z)[n(z) + f(z-u)} 

x( ; £fc^±a y + ^r + — ^) , (6b) 

- ^ + i5) + A 2 (z/ - 2 + i5; (j)') I c 2 + D 2 (u)+Q 2 (u) 



with (• • •) the angular average over 0, and 



I A(, + iM) \ 

\^Co 2 {v + i5) - A 2 (z/ + ^;0)/ 



/^ 2 (z/ + 2<5) - A 2 (z/ + i5;. 

Equations (|6]) are a set of nonlinear coupled equations for the renormalized pairing potential 
A(u + iS; 0) and the renormalized frequencies Q(y + i8) with the gap 

W + *,*) = v * v . + "-*\ (10) 

or, if the renormalization function Z(y) is introduced in the usual way as u(y + i5) = vZ(v) 
then 

A(i/ + i5; 0) = — -r . (11) 

Here, v is a real frequency and 5 is a positive infinitesimal + . To arrive at these equations 
a separable (in the angular part) model was used for the pairing potential. In the pairing 
channel it has the form g cos(20)J 2 x(u;) cos(20') with g a constant and I 2 x{^) the pairing 
spectral density. This leads to a gap proportional to cos(20) by arrangement. No other 
anisotropies are included and we note that the renormalization channel (|6b]) is isotropic 
with the same spectral density I 2 x{^) as in Eq. (|6aD but with no g value. In general a 
different form of the spectral density could come into Eqs. (|6]) but here, for simplicity, we 
have the same form but allowed for the possibility that they do not both have the same 
magnitude, i.e.: g needs not be equal to one. These equations are a minimum set and 
go beyond a BCS approach and include the inelastic scattering known to be strong in the 
cuprate superconductors. In the normal state at T near T c the inelastic scattering rate varies 
linearly in T, as it does in our model, and is of the order of a few times T c . 
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More complicated models which include the possibility of hot spotJlr§ could be intro- 
duced in our work but would not alter the main points we wish to make. In our formalism hot 
spots can be introduced by inserting an angular dependent factor in Eq. (|6b|) for the renor- 
malized frequencies. Presently, this quantity is isotropic but we could multiply the spectral 
density I 2 x{uj) (which enters this equation) by a factor which increases the scattering in the 
antinodal as compared with the nodal direction. This would complicate the numerical work 
but goes beyond what we wish to do here. In any case, recent analysis of angular resolved 
photoemission (ARPES) data on the marginal Fermi liquid (MFL) phenomenology^^ in- 
dicates that the inelastic part of the quasiparticle scattering in the cuprates may be quite 
isotropic. 

For the BCS case we get similar formulas on the real axis:^l 

A d (u,<f>) = A o cos(20) (12a) 

A» = *vrt+ / As{u) _ \ (12b) 

\ V) - A^) - A^,0) / 

u(u) = v + mt + ( . U ^ =), (12c) 

Vtf 2 M - Aft/) - Agfa / 

with impurities treated in the Born limit. 

In Eq. (|6b|) the impurity scattering rate is proportional to T + and enters only the renor- 
malization channel because we have assumed a pure <i-wave model for the gap with zero 
average over the Fermi surface. This is expected to be the case in a tetragonal system. 
The parameter c in the elastic-scattering part of Eq. (^FJ) is zero for resonant or unitary 
scattering and infinity in the Born approximation, i.e.: weak scattering limit. In this case 
the entire impurity term reduces to the form iirt + Vt(v) with c absorbed into t + . For inter- 
mediate coupling c is finite. The thermal factors appear in Eqs. (^) through the Bose and 
Fermi distribution n(z) and f(z), respectively. 

From the solutions of the generalized Eliashberg equations we can construct the 
Green's function (|3|) analytically continued to the real frequency axis u. In this formulation 
the expression for the conductivity er(T, uj) based on equations (§) and @ is lengthy; nev- 
ertheless, it is given here for completeness. For instance, the in-plane conductivity, after an 



integration over k has been performed, is given by: 



i e 2 N{0)v 2 F 
tt 2 



x ( /dz/tanh ( -^=\ —. -r —, — 

\J \2TJ E(y; 0) + E{y + Q; ■ 



+ Jdu tanh 




v + Q 



1 



2T J E*(v:6) + E*(v + Q: 



[1-N(u; 0)iV(z/ + ft;0) - 
— [1-N*(v; (j))N*(u + n; 



-P*(v; 0)P> + O; 0)] + Jdv 



1 



x 



£(z/ + fi;0) -£*(z/; 0) 


+ y dz/tanh 



z> + Q 



+ 



P*{v;<l>)P*{v + n;<f>)] 
1 



[1 + N*(u; <p)N(u + fl; 0) + 0)P(z> + fi; 0)] 
' [1 - iV*(i/; 0)JV*(z/ + fl; 0) 



£(z/ + fi;0) - E*(v;. 
with 



[1 + N*{v- 4>)N{v + fl; 0) + P*(v; 0)P(z/ + fi; . 



and 



iV(u;;0) 



(13a) 



(13b) 



(13c) 



E{u;<f>y y E(u;4>) 

In the above, (■ • •) means an average over the angle and the star refers to the complex 
conjugate. For the c-axis conductivity (■ ■ ■) is to be replaced by (■ ■ • cos 4 (20)) and the 
prefactor is different with Vp/2 to be replaced by d 2 t\. The above set of equations is 
valid for the real and imaginary part of the conductivity as a function of frequency Q. It 
contains only the paramagnetic contribution to the conductivity but this is fine since we 
will be interested in this paper mainly in the real part of the conductivity for which the 
diamagnetic contribution is zero. 

The real part of the incoherent conductivity along the c-axis is in turn given by (nor- 
malized to its normal state value cr lcTl ) illH 



Xc 



du \f(u)-f(u + 0)] 



N{u + Q)N{u) + 



V 



p(u + n)P(u) 



(14) 



s 



III. RESULTS AND DISCUSSION 



We begin our discussion by highlighting some of the important results for the real part 
of the in-plane conductivity ai(T,u) as a function of frequency at low temperatures T. We 
will want to compare our c-axis results with these in-plane results. In Fig. || we show u\{oj) 
vs. u} within BCS theory for a conventional s-wave gap (solid curve) and compare with the 
less known results for a d-wave superconducting gap function (dashed curve) .0 The solid 
curve is zero for all frequencies uj < 2Ao (twice the s-wave gap value of Ao). This is a well 
known result which has its origin in the fact that, at zero temperature, all the carriers have 
condensed into the superffuid and for there to be absorption, two quasiparticle excitations 
corresponding to the breaking of Cooper pairs out of the condensate are needed. This costs 
2Ao of energy. In the normal state (dotted curve) the equivalent excitations correspond to 
the creation of a hole and a particle (2 excitations). At u — 2A there is a sharp absorption 
edge and the conductivity jumps to a value close to its normal state Drude value (dotted line 
in Fig. |I|). In the above, the gap has been set at 24meV, a value characteristic of optimally 
doped YBCO. Also for the normal state Drude (dotted curve) the impurity scattering rate 
was t + = 0.001 meV. This very small value was employed so as to be able to resolve fine 
structures in the d-wave results which we will now describe. 

The dashed line is for a d-wave superconductor in the BCS limit. In the two dimensional 
Cu02 Brillouin zone the angle <fi gives the direction of momentum k and A^ = A cos(20) 
in our simple model. It is clear that, entering into the superconducting state reduces ai(u) 
in the region of the gap A but does not make it go to zero because even at T = some of 
the electrons have zero or very small gap values, namely those on the main diagonals of the 
Brillouin zone for = ±7r/4. It is also worth emphasizing that, above 2A the conductivity 
goes back to its normal state value faster than it does for the s-wave case. For s-wave the 
curve for ai(u) is above the normal state Drude for u > 2Ao and merges with the Drude 
form from above. For <i-wave the superconducting &i(uj) is slightly below the Drude in the 
same region but not by very much. In both cases the missing optical spectral weight on 
entering the superconducting state, of course, goes into the superfluid which provides a delta 
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function response in <J\(uj) at u = 0. This is not shown in our figure. Note also that the d- 
wave conductivity has a small structure (on the scale of the figure) at uj = A corresponding 
to a maximum in slope for a\(oS) and has another small structure at uj = 2Aq corresponding 
to a kink in the curve. These structures have their origin in the particularities of the <i-wave 
state. For example, in ci-wave, the quasiparticle density of states as a function of energy 
is linear in uj at small uj and has a logarithmic singularity at uj = A before taking on its 
normal state value at large uj. This singularity is much weaker than for the corresponding 
s-wave case for which the density of states is zero up to uj = Aq at which point it displays 
a inverse square root singularity. This square root singularity does not imply a singularity 
in <J\{oj) however. As we have seen in Fig. [I], there is instead only a sharp edge in ai(u) at 
twice the gap value. We can therefore expect much smaller structures in (Tx(u) at uj = A 
and 2A in (i-wave than in the s-wave case at 2A and this is confirmed by our numerical 
results given in Fig. [I]. 

These results are presented here mainly for comparison with full Eliashberg results based 
on the formalism given in the pervious section which includes inelastic effects. These can 
be very large in the cuprates, and are described by the charge carrier-mediating-boson 
spectral density I 2 x{uj) introduced in the previous section. Several different theoretical 
results are presented in Fig. |2|. Four different curves are shown for comparison with each 
other. They were chosen to illustrate the main features present in the in-plane conductivity 
when inelastic scattering is included. The gray solid line with a maximum before uj = 50 meV 
are BCS results as in the previous figure. It is totally different from the other curves 
and disagree strongly with experiment (solid curve). The other curves are all theoretical 
and based on the Eliashberg equations but with different models for the charge carrier- 
mediating-boson spectral density I 2 x{^). For the dash-dotted curve referred to as MMP 
with ujsf = 20meV, we employed the simple spectrum for the interaction with the spin 
fluctuations first introduced by Pines and coworkers. 00 It has the form:i!00 



where I 2 is a coupling of the carriers to the spin fluctuations and ujsf is a characteristic 




(15) 
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spin fluctuation frequency. In our work I 2 and u>sf are adjusted to get the best possible 
fit to the normal state conductivity. The same spectrum is used at all temperatures in the 
superconducting state with the value of g in the A-channel fl6a]) of the linearized Eliashberg 
equations set to get the measured value of T c = 92.4 K. There are no other parameters with 
g = 0.98 in Eqs. (j|). The resulting dash-dotted curve is to be compared with the BCS 
result. It is seen to be quite different. In particular there is no trace of the small structures 
found in the BCS (T\{oS) at u = A and 2A . In the MMP model the inelastic scattering 
at such energies is large and smears out any such small structures. However, now there is 
a significant structure at Ao plus the energy of the peak in I 2 x{oj) namely uj^y = 20meV. 
Since in a <i-wave superconductor electrons at the Fermi energy with momentum along 
the main diagonals have zero gap, the boson assisted processes start, in principle, at zero 
energy plus the boson energy {oJb) rather than twice the gap plus wg. But such processes 
make only a small contribution to the real part of the conductivity and are not seen as 
a prominent structure in this quantity. At the gap energy Ao the density of electronic 
states, however, has a logarithmic van Hove singularity and this is sufficient to produce 
the structure at A + lob described above. This is discussed in more detail in the work of 
Carbotte and coworkers&§~@ and we return to it below. It is sufficient here to state that 
a reasonable picture of the underlying I 2 x{<^) can be obtained from a second derivative of 
the corresponding optical scattering rates as a function of energy. This quantity contains 
a recognizable picture of I 2 x{^)- Here we point out that after the main Drude peak has 
largely decayed with increasing u, there is a region of low conductivity in (T\{oj) which is 
followed by a rise at higher energies. This rise corresponds to the boson assisted incoherent 
part of the electron quasiparticle spectral density. The coherent delta function which is also 
present in the quasiparticle spectral density contributes the Drude. In our work this part 
contains about 25% of the total spectral weight. 

The boson assisted region at higher energies can be used to model I 2 x{^) rather than 
taking it from some simple theory such as MMP. We can use the solid curve which shows 
the in-plane datai to get an experimentally measured I 2 x{^)- When this is done we find 
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that it consists of the 41 meV spin resonance measured in inelastic spin polarized neutron 
scattering experiments, plus a background extending up to 400 meV which we model by an 
MMP spectrum (|T5|) . The resulting I 2 x{oj) gives the dashed curve which agrees with the 
experimental data in the boson region. This good agreement provides strong evidence for 
charge coupling to the spin resonance at 41 meV which exists only in the superconducting 
state of optimally doped YBCoi below T c . The calculations were performed at T = 10 K 
which is close enough to zero. The spectral weight of the spin resonance is largest at T = 
and decreases with increasing T to vanish at T = T c . The temperature dependence of the 
resonance, measured by neutrons! is also well represented in the optical data. I There are two 
other sets of results presented in Fig. ||. The dotted curve uses I 2 x(u) in the calculations 
but also includes impurities in Born scattering with t + = 0.32 meV. This does not much 
improve the agreement between theory and experiment in the frequency region below the 
boson assisted region. The dash-double-dotted curve also includes impurities but this time 
the unitary limit is used with T + = 0.63 meV. We see that we now have an excellent fit 
to the data throughout the entire frequency range. We note the impurities hardly affect 
the boson assisted region which determines I 2 x(oj) but dominate at low u and that some 
resonant scattering must be included to get agreement in this region. 

The agreement obtained between Eliashberg theory and experiment is not limited to the 
real part of a{u). In Fig. |3| we show results for the imaginary part of c(u), more precisely 
for (jja<i((jj). The curves are labeled as in Fig. || The solid line is the data which is not well 
represented by the dash-dotted curve based on an MMP model for the spectral density. To 
get agreement it is necessary to include in an Eliashberg calculation the 41 meV resonance 
as we have done for the real part of cr(u). The pure case without impurity scattering fits the 
boson assisted region as well as the other two curves with Born scattering t + = 0.32 meV 
(dotted curve) and resonant scattering T + = 0.63 meV (dash-double-dotted curve). However, 
at lower frequencies, below 25 meV, only the curve with resonant scattering fits the data well. 
Returning to the phonon assisted region the large dip in uj<J2{uj) seen around 75-80 meV and 
reproduced by our theory is a signature of the 41 meV resonance in this quantity. The dip is 
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not present in a theory of the infrared optical conductivity based on ci-wave BCS theory.§ In 
such a theory all structures in oucr^u) vs. uo fall around or below twice the value of the gap 
amplitude Ao- The structure observed in the data at approximately 2Ao plus the resonance 
energy oj r and is thus the signature of ui T in co>o" 2 (<^)- For the MMP model (dash-dotted 
curve) there is a dip too but it falls at the wrong energy. 

We give a few more details about the second derivative technique which was used by us 
to construct the underlying spectral density I 2 x{^) from infrared data. A function W(uj) 

1 d 2 



W(u) 



00 



(16) 



2tt(1lu 2 

is constructed as a guide only. In the normal state and at low temperatures this function is 
almost exactlyS'i^i equal to the input I 2 x{^) for models based on the nearly antiferromag- 
netic Fermi liquid (NAFFL). Of course, I 2 x{&) is seen in W(uj) through electronic processes. 
But in the normal state the electronic density of states N(e) is constant and so does not 
lead to additional structures in W(u) that are not in I 2 x{uj) which would then corrupt the 
signal, if the aim is to obtain I 2 x{^) from W(u). This is no longer the case in the super- 
conducting state because of the logarithmic van Hove singularities in N(e) and these do 
indeed strongly influence the shape of W(uS) and introduce additional structures in W{uj) 
corresponding to combinations of the positions of the singularity in N(e) and the peak in 
I 2 x{^) at coy (resonance frequency) as described by Abanov et aZ.El The structures in W(uS) 
corresponding to these singularities contaminate the signal in the sense that W(u>) in the 
superconducting state is no longer equal to the input J 2 x(u;).ii In fact, only the resonance 
peak appears clearly at A + u> r and its size in W(uS) is about twice the value of I 2 x(oj) 
at that frequency. In some cases the tails in W(u) also match well the tails in I 2 x{uj). In 
the end, of course, W(uS) serves only as a guide and it is the quality of the final fit to the 
conductivity data that determines the quality of the derived I 2 x(yj). 

Nevertheless, besides giving a measure of the coupling of the charge carriers to the 
spin resonance W{u) can also be used to see the position of density of states singularities, 
as shown in Fig. £| where I 2 x{u) (gray squares) and W(uj) (solid line) derived from our 
theoretical results are compared. Also shown by vertical arrows are the positions of A + oy, 
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2A + u> r , A + 2u r , and 2A + 2u r . We note structures at each of these places and this 
information is valuable. Note that at 2A + u> r the large negative oscillation seen in W(uS) 
is mainly caused by the kink in I 2 x{^) (gray squares) at about 55meV. The density of 
electronic states effects clearly distorted the spectrum above the resonance peak and W(u) 
stops agreeing with the input I 2 x{^) in this region until about 150 meV where agreement 
is recovered. In summary, W(oj) contains some information on singularities in N(e) as 
well as on the shape and size of I 2 x{u) and, in the superconducting state, the two effects 
cannot be clearly separated. Nevertheless, W(u) remains a valuable intermediate step in 
the construction of a charge carrier-exchange boson interaction spectral density from optical 
data. 

We next turn to the c-axis dynamics and present results first for a BCS model. To 
proceed we need to specify the transverse coupling. We have presented in the theory section 
two possible models. The first one is coherent tunneling with the matrix element £j_(k) 
conserving in-plane momentum and probably equal to t± cos 2 (20) where is the angle of k 
in the two dimensional CuC>2 Brillouin zone. The other model is incoherent tunneling for 
which momentum is not conserved. To be definite we will use |Vi/Vo| = 1 ( see Eq. (|5|)) in 
the impurity potential. Other values of |Vi/Vo| have been studied by Hirschfeld et and 
by others!!! to which the reader is referred for more details. 

In Fig. |5| we show BCS results for the real part of the c-axis conductivity (Ti c {oS) vs. uj and 
compare with the in-plane conductivity cti(uj) (dotted curve). In this case the gap is 24meV 
and the in-plane impurity scattering rate is t + = 0.1 meV. We see that the structure at A 
and 2A in the dotted curve for <Ji(u) have been smeared out somewhat more when compared 
with the equivalent results shown in Fig. |l[ This is due to the larger impurity content. This 
curve is for comparison with our c-axis results which we now describe. The solid curve is for 
the coherent tunneling case with ij_(k) = tj_cos 2 (20). This hopping probability eliminates 
the nodal quasiparticles along (it, tt) which do not participate in the out-of-plane dynamics. 
This gets rid of much of the remaining Drude like contribution at very low u which is still 
clearly present in the dotted curve although it is substantially reduced by superconductivity 
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as compared to the normal state Drude (see Fig. [I]). The solid curve is small at small uj 
and peaks just below uj ~ 2A where the dotted curve for the in-plane case has a small 
structure. At higher energies, there is little difference between the in-plane and out-of-plane 
<Ji{u) although the magnitude of these two quantities is of course very different. If instead 
of using £j_(k) we had used a constant t± for the c-axis transport, the out-of-plane <J\ c {uj) 
would mirror the in-plane case (?i(uj). This would also hold in the more complex Eliashberg 
calculations. Only the magnitude is different between in-plane and out-of-plane in this case 
because of differences in over all multiplicative factors in front of the expression for the 
conductivity. The final curve in Fig. |^, dash-dotted, is for the incoherent case. We see 
that it too is near zero at small uj although it rises out of zero more rapidly than does the 
solid curve. It shows no structure whatever at the gap or at twice its value. The main 
rise is accomplished within the region uj ~ 2A . At high uj it saturates to a constant value 
of one. This is because we have normalized our results to the normal state conductivity 
and the curve becomes very flat. The corresponding normal state conductivity would be a 
horizontal line at this saturated value, constant for all uj. This behavior bares no relation 
to the in-plane coherent result. 

The results for ai c (uj) presented so far are for comparison with those based on solutions for 
the Eliashberg equations given in the previous section and which properly include inelastic 
scattering through the spectral density I 2 x{oj). Before presenting our c-axis results in this 
case we stress again that the boson exchange kernel I 2 x{u) is an in-plane quantity and is 
taken from our discussion of the in-plane conductivity. It is not fitted to any c-axis data. It is 
to be used unchanged to calculate the out-of-plane conductivity assuming coherent hopping 
with £j_(k) = t± cos 2 (0). The solid curve in Fig. |6] are the in-plane Eliashberg results which 
are included for comparison with the dashed curve which is for the c-axis. In the boson 
assisted region, which would not exist in a BCS theory, both curves have a remarkably 
similar behavior. At very low frequencies, a region which comes mainly from the coherent 
delta function part of the carrier spectral density, and which is the only part included in BCS, 
we note a narrow Drude-like peak in the solid curve. This part is suppressed in the c-direction 
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(dashed curve) because the contribution from the nodal quasiparticles are effectively left out 
by the t± cos 2 (0) weighting term. Also, shown for comparison are our previous BCS results 
for coherent hopping (dotted curve). These results show no resemblance to our Eliashberg 
results and also do not agree with experiment. What determines the main rise in the region 
beyond the Drude part of the conductivity in cr lc (u;) are the boson assisted processes and 
this rise does not signal the value of the gap or twice the gap for that matter but rather a 
combination of A and the resonance energy u r . 

In Fig. ^ we compare the data from Homes et on the same graph for in-plane 
(dotted) and out-of-plane (solid) conductivity ai(o>). It is clear that in the c-direction, 
the nodal quasiparticle seen in the dotted curve are strongly suppressed. This favors the 
t j_ cos 2 (2(f)) matrix element for the c-axis dynamics as we have just seen. Further, in the 
boson assisted region the two curves show almost perfect agreement with each other, which 
again favors the t± cos 2 (20) coupling as illustrated in the theoretical curves of Fig. || One 
difference is that the main rise, indicating the onset of the boson assisted incoherent (in- 
plane) processes, appears to have shifted slightly toward lower frequencies in the c-axis data 
as opposed to a shift to slightly higher frequencies in our theory. It should be remembered, 
however, that in the raw c-axis data, large structures appear in the conductivity due to 
direct phonon absorption and these need to be subtracted out, before data for the electronic 
background of Fig. [7| can be obtained. In view of this, it is not clear to us how seriously we 
should take the relatively small disagreements that we have just described between theory 
and experiment. 

With the above reservation kept in mind we show in our last Fig. |8| a comparison of various 
theoretical results with experimental c-axis conductivity (black solid line). There are five 
additional curves. The black ones are obtained from an Eliashberg calculation based on the 
MMP model for I 2 x{^) with impurities t + = 0.32 meV included to simulate the fact that the 
samples used are not perfect i.e. are not completely pure, but this parameter does not play a 
critical role in our discussion. Incoherent c-axis coupling is assumed with |Vi/Vo| = 1 (black 
dotted). It is clear that this curve does not agree with the data and that the coupling along 
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c cannot be dominated by incoherent hopping between planes. This is also in agreement 
with the results of a theoretical study by Dahm et a/.@ who also observed better agreement 
for coherent c-axis conductivity in the overdoped regime. On the other hand the fit to the 
black dashed line is good in comparison. It uses the same MMP model but with coherent 
coupling of the form t±(k) = t^cos 2 (20). This fit may already be judged satisfactorily 
but it should be remembered that if we had used the model of 7 2 x(u;) with the 41meV 
peak included instead of MMP, the agreement would have deteriorated. This is troubling 
since one would expect that coupling to the 41 meV spin resonance would be stronger in 
the c-direction data than it is in the in-plane data. This is because the c-axis emphasizes 
the hot spots around the antinodal directions which connect best to (tt, tt) in the magnetic 
susceptibility. This is the position in momentum space where this spin resonance is seen to 
be located in optimally doped YBCO. On the other hand, recent ARPES data&B which 
fit well the MFL (marginal Fermi liquid) phenomenology show little in-plane anisotropy for 
scattering around the Fermi surface and this is consistent with the findings here. 

The dash-dotted curve in Fig. ^ illustrates the fit to the data that can be achieved with a 
dominant coherent piece and subdominant incoherent contribution. It is not clear to us that 
such a close fit is significant given the uncertainties in the data and the lack of uniqueness in 
the fitting procedure. It does, however, illustrate the fact that a small amount of incoherent 
c-axis hopping cannot be completely ruled out from consideration of the infrared data and 
that this data can be understood quite well within Eliashberg theory. The last two curves 
(solid gray and dotted gray) are based on BCS d-w&ve theory and are reproduced here to 
illustrate the fact that such a theory is unable to explain the c-axis data. The solid gray 
curve is with t±(k) = tj_cos 2 (20) and the dotted gray one for incoherent c-axis transport. 
Compared with our Eliashberg results the agreement with the data is poor. 

IV. CONCLUSION 

We have considered the c-axis charge response as revealed in the real part of the frequency 
dependent optical conductivity cr lc (c<;). Results for a pure BCS model, which includes only 
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elastic impurity scattering, show no agreement with the data in- or out-of-plane. A gen- 
eralized Eliashberg approach based on a spin fluctuation mechanism leads to much better 
agreement. 

In this approach the coherent delta function like part of the electron spectral density 
which is sharply peaked at the quasiparticle energy leads to a Drude type response in 
the normal state and a BCS type response in the superconducting state. The dominant 
incoherent background of the electron spectral density provided by the coupling to the 
spin fluctuations, however, leads to an additional boson assisted region for the conductivity 
which is not strongly changed by the onset of superconductivity except that it is a shift 
by roughly Ao- This part of the conductivity, not included in BCS, contains the largest 
part of the optical spectral weight (of order 75%) and can be used to get information on 
the underlying inelastic processes that presumably cause superconductivity. The in-plane 
conductivity reveals coupling to the 41 meV resonant peak seen in spin polarized neutron 
scattering as well as to the spin fluctuation background extending to high energies. For 
coherent c-axis hopping the boson region in ai c (u>) has the same form as it does in the 
in-plane case and this fact is largely born out in the experimental data. This would not be 
the case if the c-axis coupling were incoherent. Assuming pure incoherent coupling gives no 
agreement with the data although a best fit is obtained with a small subdominant incoherent 
part in addition to a dominant coherent part. An important conclusion coming out of our 
analysis is that the c-axis conductivity data gives independent confirmation for the form of 
the charge carrier boson spectral density obtained solely from in-plane infrared optical data. 

While the infrared conductivity at higher energies is dominated by the inelastic processes 
and not much affected by a modest amount of elastic impurity scattering, the low energy 
part is much more sensitive to impurities and depends less directly on inelastic scattering. 
In this sense, this frequency region can be partly understood within BCS theory at least 
at low temperatures. One needs to remember, however, that the spectral weight under the 
Drude-like curve is only a fraction of the entire spectral weight and comes only from the 
delta function-like part of the electron spectral density. Also, even in this energy region, 
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there are other profound modifications introduced by the inelastic scattering off the spin 
fluctuation spectrum. For example all traces of the small structures expected in a\{uS) at 
the gap and twice the gap values are smeared out because, at such higher frequencies, the 
inelastic scattering rate is already large. Structures do appear in the data, however, at the 
energy of the gap and of twice the gap plus once or twice the energy of the spin resonance 
peak which is seen by inelastic neutron scattering at 41meV in optimally doped YBCO. 

Finally, we point out that the spectral density for the excitation spectrum that causes 
the superconductivity in our approach which consists of a peak at 41 meV plus a long 
nearly constant background extending to very large energies while consistent with NAFFL 
of Pines and coworkers, 00 could also arise in other microscopic mechanisms. Recently, the 
marginal Fermi liquid model, proposed early on, in the development of our understanding 
of the cuprates, and found to reproduce very well many of the observed anomalous normal 
state properties, has received new attention because it also fits well the ARPES data. The 
excitation spectrum of the MFL is quite similar to the one of NAFFL and has many of 
the features of our empirically determined spectrum provided its normal state version is 
modified by the 41 meV resonance on entering the superconducting state. In that sense the 
MFL is equally consistent with optical data. Both, MFL and NAFFL, are phenomenological 
models having similar excitation spectra associated with the pairing and so cannot easily be 
distinguished from optics. 
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FIGURES 

FIG. 1. Comparison of the real part of the in-plane conductivity (J\(u) of a d-wave supercon- 
ductor (dashed line) and an s-wave superconductor (solid line). Each have a gap of Ao = 24meV. 
For the s-wave there is no absorption till 2Aq while for the d-wave with A^ = Ao cos (20) with cj) 
an angle in the two dimensional Brillouin zone, there is a reduction as compared to the normal 
state Drude (dotted curve). The slope of (Ti(uS) has a maximum at Ao and cri(u) has an additional 
small structure at 2Ao- 

FIG. 2. Comparison of the real part of the in-plane conductivity (J\(u) vs. uj for various models. 
The grayed solid line with a peak before 50meV is BCS. The dash-dotted line is an Eliashberg 
calculation with an MMP spectral density peaked at wgp = 20 meV. The dashed line is the same but 
with I 2 x(u) used instead of the MMP model. As described in the text this electron-boson spectral 
density I 2 x(u) has been determined through a consideration of the in-plane optical data. The 
dotted (Born) and dash-double-dotted (unitary scattering) curves include impurities in addition 
to the I 2 x(w) model for inelastic scattering. The solid line is the data of Homes et aZ.HH 

FIG. 3. The imaginary part of the conductivity lu(J2(uj) vs. u for the various models described 
m Fig. g. The solid curve is the data. The dash-dotted curve is the result of an Eliashberg 
calculation with an MMP model while the other curves are based on the model I 2 x(u) which 
includes the 41 meV resonance. These three curves are for the pure case (only inelastic scattering, 
dashed line), and with some additional elastic impurity scattering in Born (dotted) and unitary 
(dash-double-dotted) limit with t + = 0.32 meV and r + = 0.63 meV respectively. 

FIG. 4. Second derivative W(uj) compared with input boson spectral density I 2 x(w). The 
41 meV peak in I 2 \(uj) (gray squares) is clearly seen in W(lo) (solid line) as are the tails at higher 
energies. In the energy region between 75 and 150 meV the van Hove singularities in the electronic 
density of states show up added on to uj t and distort the correspondence between W(u>) and I 2 x{u). 
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FIG. 5. Comparison of the real part of the conductivity a\(uj) vs. lo for the in-plane case 
(dotted curve) and c-axis. The solid line is for coherent c-axis hopping with t((p) = ij_ cos 2 (</>) 
with <f) an angle in the two dimensional CuO"2 Brillouin zone. The dash-dotted is for incoherent 
tunneling with |Vi/Vo| = 1. All curves are in the BCS model. 

FIG. 6. Comparison between (solid) in-plane and out of plane (dashed) real part of the c-axis 
conductivity ai c {uo) vs. lo in an Eliashberg model with our model carrier-boson spectral density 
I 2 x(lo) which includes the 41 meV spin resonance. The dotted curve is ai c (uj) for a BCS d-wave 
model with the same gap value as in the Eliashberg work and is included for comparison. 

FIG. 7. Comparison between in-plane (dotted) and out of plane (solid) for the real part of 
oi{u) vs. u. The data is from Homes et aL0 

FIG. 8. Comparison with the data of Homes et a/.0 for the c-axis conductivity (black solid 
curve). The theoretical curves were obtained in a BCS theory, solid gray (coherent), dotted gray 
(incoherent) and the others in Eliashberg theory with MMP model and impurities t + = 0.311 meV. 
The black dotted curve is for incoherent c-axis with |Vi/Vo| = 1, the dashed for coherent c-axis 
with t±((j)) = t± cos 2 (2^) with <fi an angle in the two dimensional Cu02 Brillouin zone, and the 
dash dotted is a fit to the data provided by a mixture of coherent and incoherent. We stress that 
this last fit is for illustrative purposes only, and is not unique. 
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